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Abstract We investigate renormalization group limit cycles within the similarity renormalization 
group (SRG) and discuss their signatures in the evolved interaction. A quantitative method to detect 
limit cycles in the interaction and to extract their period is proposed. Several SRG generators are 
compared regarding their suitability for this purpose. As a test case, we consider the limit cycle of the 
inverse square potential. 
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1 Introduction 

Few-body systems with resonant interactions have universal properties independent of the details of 
the interaction at short distances [T]. If the scattering length a is much larger than the range of 
the interaction ro, the Efimov effect can occur [2]: there is a geometric spectrum of three-body bound 
states with an accumulation of states at threshold. For identical bosons in the unitary limit, the Efimov 
spectrum satisfies the geometric scaling relation 

/ \ 2 

where 5o = 1.00624... is a transcendental number. In general, the value of sq depends on the sym¬ 
metries of the three-body system, the number of interacting pairs, and the masses of the particles 
involved, but the form of the spectrum 0 is universal. The invariance of the spectrum under discrete 
scale transformations with the preferred scaling factor squared (22.7)^ can be understood as 

the consequence of a three-body interaction H governed by a renormalization group (RG) limit cycle 
OSIIS]. The coupling constant H{A) is then a periodic function of In(yl) with period tt/sq where A is 
an ultraviolet cutoff in momentum space used to regulate the theory at short distances. More generally, 
the limit cycle will be manifest in observables through their log-periodic dependence on the scattering 
length or other control parameters of the system [T] . The period of this log-periodic behavior is again 
determined by the discrete scaling factor which is the key quantity determining the properties of the 
limit cycle. 

The universal Efimov spectrum and related higher-body bound states have been observed in ul¬ 
tracold atomic gases in a variety of experiments with different atom species miziiH]. In heteronuclear 
mixtures the scaling factors can be significantly smaller than for ideal bosons. For example, in a mix¬ 
ture of ^Li and ^^^Cs atoms, the ratio of subsequent ^^^Cs-^^^Cs-^Li bound states is predicted to be 
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(4.88)^. Such heteronuclear Efimov states have recently been observed in experiment and the predic¬ 
tion of the scaling factor was confirmed muni. Finally, we note that the Efimov effect also provides 
a universal binding mechanism for S-wave states in weakly-bound nuclei. However, an observation of 
the discrete scaling relation between different states in one nucleus has proven elusive to date since the 
known S-wave halo nuclei have no excited states. (See Ref. [11] for a review.) 

Here, we investigate the manifestation of RG limit cycles in the similarity renormalization group 
(SRG). The SRG was developed independently by Wegner [12] and by Glazek and Wilson [13]. In 
recent years, it has become a standard tool to soften nuclear interactions for improved convergence in 
many-body calculations [ims]. The SRG generates a continuous series of unitary transformations of 
the Hamiltonian from the evolution with a flow parameter s. As the RG flow evolves, the interaction 
is typically softened at the expense of introducing induced many-body forces. Ideally, the evolution is 
carried out to a value s such that the interaction is softend enough to achieve convergence in many-body 
calculations without generating large contributions from induced many-body forces. Recent advances in 
SRG technology have, e.g., allowed to extend ab initio calculations to P-shell nuclei [16] or perturbative 
neutron matter calculations with consistent three-body forces HU. The success of SRG methods in 
nuclear physics motivates our study of the manifestation of RG limit cycles within the SRG. 

The emphasis of this work is on formulating criteria for detecting a limit cycle in the evolved 
interaction^ In particular, we study the numerical extraction of the discrete scaling factor. A long 
term goal of our work is the derivation of the limit cycle in the pionless EFT for the nuclear three- 
body system. While it is generally assumed that the pionless EFT in which the Efimov effect and 
corresponding limit cycle are manifest appears as the low-energy limit of a more fundamental chiral 
EFT with explicit pions, no explicit derivation has ever been given. Our work provides a first step 
towards this goal. As a test case, we investigate the attractive inverse square potential which is known 
to have an exact limit cycle beyond a critical coupling strength. For convenience, we will use natural 
units with h = m = 1 in the following. 


2 SRG Basics 


We start with a brief review of the basic properties of the SRG. A more detailed discussion can be 
found in Refs. [muH]. 

The SRG generates a continuous series of unitary transformations on the Hamiltonian governed by 
a flow parameter s: 

H{s) =U{s)H{s = 0)U{s)^ =T + V{s) , (2) 

where U{s) is a unitary operator. Often the kinetic energy is left unchanged per definition, such that 
the SRG generates an evolution of the interaction potential V. Defining the anti-Hermitian operator 
r]{s) = {-^U{s))U{sy = [G{s), H{s)], the evolution equation for the Hamiltonian can be written as 

^ = [[G{s),H{s)],H{s)] , (3) 

where the Hermitian operator G{s) is called the generator of the SRG transformation. The generator 
G is often taken independent of s. A frequent choice for G is the kinetic energy T. In this case, the 
differential equation can be written as 


A = 2TVT - VTT - TTV + TVV + VVT - 2VTV . 

as 


( 4 ) 


If we consider a two-body system with identical particles, we can write the flow equation for the 
two-body interaction V 2 in the space of relative momenta as 


■^{p\V2\q) = - q^f (p|V'2 


k) + r 

Jo 


d^k 


ip^ + q'^ -2k‘^){p\V2\k){k\V2\q) 


( 5 ) 


^ We note that a limit cycle will also be manifest in physical observables through universal scaling relations 


m- 
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The evolution of partial waves decouple in the two-body system, thus Eq. (§ holds for every partial 
wave. In this representation one can recognize a major characteristic of the SRG with the T generator. 
Clearly, the SRG transformation has a fixed point if H{s) commutes with T, i.e. H{s) is diagonal. 

For weak interactions, the second term in Eq. (§ can be neglected and the solution is simply 

{p\V2{s)\q) = (p|V2(0)|g)exp(-s(p2 - q^f) . (6) 

All non-diagonal matrix elements approach zero during the evolution and thus a decoupling between 
low-energy and high-energy matrix elements is achieved. Although the second term in Eq. (§ typically 
can not be neglected, the suppression of non-diagonal matrix elements is a dominant part of the full 
SRG evolution using the T generator. The approximate expression ^ also suggests to interpret 

A = (7) 


as an effective momentum cutoff. Thus the evolution starts at 5 = 0 corresponding to A = oc and as s 
increases, the effective momentum cutoff A is lowered. 

In this work, we also use two alternative generators, which are functions of the kinetic energy. They 
were introduced by Li, Anderson, and Furnstahl with the aim to obtain a more efficient evolution 
which is crucial for identifying limit cycles m- We refer to them as exponential generator, Ge, and 
inverse generator, Gf. 

Ge = -o-^exp(-r/(T^), 

^ 1 + T/(t2 ’ 


where a is an arbitrary parameter with dimensions of momentum. Both generators have a power series 
expansion in T. For small momenta q a, they approach T up to a constant which cancels out 
in 7^. So depending on the parameter a, there is a separation into a low-energy region, where the 
two generators behave like the T generator and a high-energy region, where the SRG evolution is 
suppressed. A detailed discussion of these generators can be found in Ref. [19]. 

Since a has dimensions of momentum, the translation of s to an effective momentum cutoff is more 
subtle than for the T generator. The solutions of Eq. in the weak interaction limit are 


(pV 2 (s)k) = (p|V' 2 ( 0 )|g)exp[-S(T 2 (g 2 -p2)(e - e , 


for Ge , 


( 9 ) 


and 


{p\V2{s)\q) = (p|V2(0)|g)exp 


-saHq^-p^) 


1 + 1 + g^/c 


for Gj . 


( 10 ) 


Thus the effective momentum cutoff is A ~ a The constant a is irrelevant for the parametric 

dependence on s. Setting it to one, we define the effective momentum cutoff 


An, — S 


- 1/2 


( 11 ) 


for the exponential and inverse generators from Eq. Q. This will have important consequences when 
extracting limit cycle periods from the evolved interaction below. 


3 Renormalization of the 1/R^-Potential 

In the following, we discuss the quantum mechanical 1/R‘^ potential as a test case. This is a singular 
potential which displays an exact limit cycle. We start by reviewing the renormalization of the 1/R^ 
potential in an effective field theory framework. Here, the limit cycle becomes manifest in the behavior 
of a counter term. We follow the discussion in Ref. [20] where further details can found. In the next 
section, we investigate the 1/R^ potential in the SRG framework and provide general criteria for 
isolating limit cycle behavior in the interaction. 

The 1/R^ potential can be written as 

V'(fl) = ^ 


(12) 
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with R := |R| and c a coupling constant. For subcritical couplings c > — the potential is well 
behaved and leads to a unique solution of the Schrodinger equation. For critical and supercritical 
values c < — |, however, the potential is singular and displays a limit cycle. In this case, it is useful to 

define a parameter u that characterizes the period of the limit cycle via u := c — 

The momentum-space representation of the potential can be defined via a Fourier transform in D 
dimensions, 

27r^c 


V{Q) = lim 

D^3 


/ 


^DRgiQ.R 


ViR) = 


Q 


(13) 


where Q is the momentum transfer [20] . 

In the following, we consider only S-waves. For the momentum space matrix elements of the S-wave 
projected potential, we get 

'0{p-q) 0{q-p)' 


V (p, q) = 27r^c 


(14) 


p q 

where q (p) are the incoming (outgoing) momenta. The physical observables can be obtained from the 
Lippmann-Schwinger equation 

... 

:V{p,q)tE{q,k) , (15) 


tE{p,k) = V{p,k) + ^ 


- f 

Jo 


dqq^ 


■ le 


E-q^ 

where E = the total energy and the scattering phase shifts are given by 

47r 

k cot S = ik - - ——-. 

tE{k,k)\^^j^2 

The bound states are given by the solutions of the corresponding homogeneous equation. As discussed 
in [20], Eq. 
regulate Eq. 


(16) 


(| l5|) has no unique solution for yl —> oo if c < — j and requires renormalization. We 
with a sharp momentum cutoff A and absorb the cutoff dependence by introducing 


a momentum independent counterterm SV (yl). 

V{p, q) => V{p, q) + SV{A) = 27r^c 




p q A 

Demanding invariance of the zero-energy solution under changes of A, one finds 


H{A)^^ 


H{A) = 


1 — 2utQon{yhi{A/A^)) 

r 


= 1 - 4z/^ ln(A/A*) + . 


(17) 


(18) 


■ 2z/tan(z/ln(A/A^)) 

where A^ is a low-energy constant. Including this counterterm in Eq. ( [T^ keeps all low-energy ob¬ 
servables fixed when A is varied. One can immediately see that the counterterm E[(A) displays a limit 
cycle with a preferred scaling factor exp(7r/z/) since tan is a periodic function with period tt. If the 
cutoff A is changed by multiples of exp(7r/z/), the counterterm returns to the same value. 

The bound state spectrum satisfies a geometrical scaling relation analog to the Efimov case, 

-E- _= fl9) 

p;(n+i) ^ 

and presents an ideal test case for the application of SRG methods to limit cycles. In the following 
section, we will investigate the limit cycle in the inverse square potential using the SRG framework. 


4 XjR? Potential and SRG 

In this section, we consider the 1/R^ potential in the SRG framework. Since the SRG is a unitary 
transformation, all observables stay constant during the evolution by definition. So in contrast to 
the explicit construction of the counterterm in the effective field theory treatment of the previous 
subsection, we need to extract a signal for the limit cycle from the evolved interaction. In order to 
define such a signal, we investigate the SRG evolution of the 1/R^ potential for critical and subcritical 
couplings and different generators. 

Before we proceed, we specify our units. One free length scale Iq is present in our framework. 
Therefore, all dimensionful quantities are given in units of Iq. 
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Fig. 1 Evolution of the S-wave potential for the T generator with z/ = 9 and initial cutoff yl = 20 

V{k,k',X) is shown for A ~ 21.09, 13.67, 8.86, 5.74, 3.72 and 2.41 in units of [l//o] from top left to bottom 
right. 


4.1 Qualitative Features 


First, we consider the standard T generator for the SRG transformation. In Fig. 
evolution of the potential for = 9 and an initial cutoff in Eq. (15) of yl = 20 /, 
regulator is required in order to insure that Eq. (15) has a unique solution. 


7-1 


we show the 
Introducing this 
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k[i/y k-[i/y 



Fig. 2 S-wave 1/i?^ potential evolved with the T generator to A = 1 /q ^ with initial yl = 20 /q ^ • The potential 
strengths are c = —0.125, c = —0.25, c = —0.25 — and c = —0.25 — 2^ in the order from top left to bottom 
right. 


Changing the value of A corresponds to changing the short-distance behavior of the starting in¬ 
teraction. We leave A fixed and investigate the dependence of the interaction on the SRG cutoff A 
that is directly related to the flow parameter s. Striking is the appearance of separated regions in the 
potential with positive and negative sign. To make them clearly visible, we choose a rather small max¬ 
imum value for the coloring of the potential. As the evolution progresses, these regions are constantly 
emerging and vanishing while the total number of regions increases. At the beginning of the evolution, 
two positive valued regions appear. During the further progress more and more positive and negative 
regions emerge. The size of these structures also decreases, which is related to the general suppression 
of off-diagonal matrix elements in the SRG evolution for the T generator. In the last picture, one can 
clearly see the large number of small regions. We also note that this behavior occurs on a logarithmic 
scale of the flow parameter s. 

To confirm that the appearance of the oscillatory behavior is indeed related to the limit cycle, we 
have evolved 1/R^ potentials with subcritical coupling c> where no limit cycle occurs and critical 
couplings c < —\ where it is present. In Fig. four potentials are depicted, which were all evolved 
to A = 1 All parameters of the potentials except for the coupling constant c are kept constant. 

For c < — ^, the scaling factor is given by exp(7r/z/) with v = c — Thus, if c approaches the 

critical value — \ the scaling factor diverges. In Fig.|^ the evolved potentials beneath and at the critical 
value do not exhibit the oscillatory behavior. Only the effective diagonalization of V from the SRG 
transformation is clearly visible. For = 1, two separated regions with opposite signs are observable 
up to this point in the evolution and for z/ = 2 several structures are already visible. This observation 
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Fig. 3 Evolved potential with parameters z/ = 9, /I = 20 and a = 2 for Ge and Gi. Left 

panel: Evolution with the exponential generator Ge to Aa ~ 8.35 Right panel: Evolution with the inverse 
generator Gi to Aa ~ 10.75 


clearly supports our conjecture that the appearance of the oscillatory features is related to the limit 
cycle. 

Next, we consider the inverse and exponential generators. We expect a similar qualitative signature 
of the limit cycle. However, the alternative generators contain a free dimensionful parameter a which 
divides the potential into two different regions. For small momenta compared to a, the exponential and 
inverse generators reduce to the T generator. For large momenta, the generators approach zero and the 
evolution is suppressed. As an example, we have plotted the evolved potential for = 9, A = 20 
cr = 2 in Fig. for both the inverse and the exponential generator. 

Notable is the fact that the oscillatory features become compressed in a rather small area in the 
k — k' plane, whose size depends on a. We remark that the qualitative behavior of the exponential and 
inverse generators are very similar. Thus, we will not distinguish their traits here. As for the standard 
T generator, the oscillatory structures only appear if the coupling constant is supercritical. 


4.2 Discrete Scaling Factor 


Next we focus on the question of how to extract the discrete scaling factor as the main characteristic 
of the limit cycle from the evolved potential. To observe a log-periodic signal, we investigated several 
different strategies, which we will discuss in the following. 

First, we have investigated the possibility to isolate the oscillatory feature by projecting on the 
momentum independent part of the evolved potential. Assuming that the short-range part of the 
evolved potential can be expanded as 

V{k, k', A) = Co(A) + + fc'2) + ... , (20) 

we have investigated the possibility to extract the preferred scaling factor from the A dependence of 
Cq. In particular, we considered the quantities 


Zi(A) 


(fk(fk'V{k,k',X) 


and 



d^kV{k,k,\) . 


( 21 ) 


/i(A) is the projection of the evolved potential whereas /2(A) is the projection of the diagonal part of 
the evolved potential. In both cases, we were not able to detect any clear signals of the limit cycle. 

Second, we examined the diagonal elements of the evolved potential s) in dependence of the 

flow parameter. This is motivated by the diagonalizing factor exp(—from Eq. On the 
diagonal of the potential matrix incoming and outgoing momenta are identical, so that the exponential 
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Fig. 4 (a) Diagonal element V{p,p,X) in dependence of X for p ^ 0.84 Iq (b) V'(p,p,X) from Eq. (22) in 
dependence of A for p 0.84 . The parameters of the potential are z/ = 9 and T = 20 The evolution 

was carried out with the standard T generator. 



Fig. 5 Evolution of the diagonal elements V{p,p, Xa) with z/ = 11, 4 = 30 /q ^ and a = 0.05 Iq ^. (a) Applied 
exponential generator with p ^ 1.19a (b) Applied inverse generator evolved with p ^ 1.55cr. 


function is one. Hence, the diagonal elements are the only ones, which do not approach zero during the 
evolution. So, we expect a log-periodic signal to be most prominent on the diagonal. A similar strategy 
was followed by Glazek [21] in the analysis of a discrete model displaying a limit cycle. 

T generator. We start with the standard T generator. In Fig. ia) a typical diagonal element is 
depicted in dependence of the flow parameter. The diagonal elements show some irregular oscillations 
but a clear signature of the limit cycle period can not be extracted. This is also the case if we subtract 
the initial potential from the evolved potential and multiply by A in order to isolate the SRG analog 
of the counterterm H from Eq. 0 above: 

V\p,p,X) = {V{p,p,X)-V{p,p,\ = oo)) ■ X, where A = . (22) 


In Fig. [^b), we show W(p,p, A) from Eq. ( [^ for the same diagonal matrix element as in Fig. [^a). 
Again a clear signature of the limit cycle period could not be extracted. This was also the case if 
different powers of A were used in Eq. (22). 

Alternative generators. Using the standard T generator only very few irregular oscillations can 
be observed. Therefore, we try the same strategy with the exponential and inverse generators that 
allow for a further evolution in s. Here, a completely different behavior is found if momenta of order 
a are considered. In Fig. we plot a diagonal element with momentum p close to the parameter a in 
dependence of the flow parameter for both generators. One can now clearly see regular oscillations for 
both generators. The graphs look like a log-periodic function multiplieded with another slowly-varying 
function. 

An example with explicit values is given in Table The extracted distances between the maxima 
and minima are constant to about 5%, except for the first few oscillations at large values of A^ which is 
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oscillation 

u = 

maxima 

11 

minima 

u - 

maxima 

= 5 

minima 

1 

2.94 

3.08 

3.26 

3.36 

2 

1.76 

1.80 

2.11 

2.11 

3 

1.52 

1.54 

1.91 

1.91 

4 

1.42 

1.43 

1.87 

1.85 

5 

1.38 

1.39 

1.83 

1.81 

6 

1.36 

1.36 

1.85 

1.80 

7 

1.35 

1.34 

1.84 

1.79 

8 

1.31 

1.32 


1.78 

9 

1.33 

1.33 



10 

1.34 

1.33 



11 

1.30 

1.31 



12 

1.33 

1.31 



13 

1.34 

1.34 



14 

1.38 

1.35 



15 

1.36 

1.36 



16 

1.33 

1.31 



17 


1.34 




Table 1 Examples of extracted ratios Aa for the exponential generator, where the Aa^ are the flow 

parameter values of the maxima (minima) of E(p,p, Aa) with p ^ 1.19(7 for z/ = 11 and z/ = 5. The parameter 
values are A = 30 and a = 0.05 The exact scaling factors are exp(7r/ll) ~ 1.33 and exp(7r/5) ~ 1.87. 



probably caused by finite cutoff effects. The period depends on the strength of the initial potential u. 
Larger values of u result in smaller periods. We find that the extracted periods are in good agreement 
with the exact formula exp(7r/z/). The agreement is better for larger values of z/, where more oscillations 
can be seen and the period can be determined more accurately. 

We will now elaborate on the appearance of the oscillations on the diagonal. To this effect, the 
diagonal elements V{p^p,Xa) are displayed in Fig. as a function of A^ for four different momenta. 
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The clearest oscillations can be extracted in the region p ^ a. For larger momenta, the amplitudes 
of the oscillations become smaller, cf. Fig. |^c) and (d). If p is even further increased, the oscillations 
disappear. Figure]^ (a) demonstrates that choosing momenta smaller than a leads to fewer oscillations 
on the diagonal. If p is even further decreased, the number of oscillations with large amplitude is reduced 
and the graph resembles the diagonals of potentials evolved with the T generator. This behavior is 
expected since the inverse and exponential generators reduce to the standard T generator for small 
momenta p. We briefly return to Fig. [^b). For smaller A^, the amplitude of the oscillation decreases 
until no oscillation is visible anymore. The number of observable oscillations strongly depends on a. 
Typically, they can be seen until A roughly equals a. 

In summary, we find a clear signal of the limit cycle in the A^ dependence of V(p,p, Xa) for p ^ cr. 
In this region, the evolution of the alternative generators is distinctively different from the standard 
T generator but not heavily suppressed as for large p a. In practice, many oscillatory features are 
compressed in the region p ^ a and thus can be distinguished from other effects of the evolution. They 
can be used to extract the preferred scaling factor numerically. 


5 Summary and Outlook 

In this paper, we have investigated the renormalization group limit cycle of the quantum-mechanical 
1/R^ potential within the similarity renormalization group framework. We showed that the period 
of the limit cycle can be extracted from the A^ dependence of the diagonal elements of the evolved 
interaction A^) if the alternative generators in Eq. <§ with a dimensionful parameter cr are 

used. In the region p ^ a sufficiently many oscillatory features are present, such that a numerical 
extraction of the preferred scaling factor is possible. Here, we have determined the period simply by 
taking averages over the positions of minima and maxima in E(p,p, A^). More elaborate statistical 
analysis schemes using Bayesian statistics can help to improve the extraction of the period [22]. 

In the future, it will be interesting to apply our technique to the nuclear three-body system in order 
to derive the limit cycle in the pionless EFT from a more fundamental chiral interaction with explicit 
pions. Work in this direction is in progress. 
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